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Abstract 

We show that lattice regularization of noncommutative field theories can be used 
to study non-perturbative vacuum phases. Specifically we provide evidence for the 
existence of a striped phase in two-dimensional noncommutative scalar field theory. 
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Introduction 



The mixing of UV and IR physics in noncommutative field theories is a intrigu- 
ing phenomenon. One of the more spectacular manifestations of this mixing is 
the appearance of a so-called striped phase which breaks translational invariance 
in noncommutative 4 theory. The prediction of this phase was based on a self- 
consistent Hartree treatment of one-loop diagrams [1] or a one-loop renormalization 
group approach [2]. The purpose of the following note is to test the predictions 
using the non-perturbative framework of lattice computer simulations and at the 
same time test the viability of studying noncommutative field theories by numeri- 
cal means. Here we consider the simplest possible such theory: a two-dimensional 
4 theory. Two-dimensional chiral and Yang-Mills theories have also been studied 
recently [3, 4]. 



The model 

Our starting point is noncommutative scalar field theory with the action 
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where 0* 4 = 0*0*0*0 and where the star product is defined by 
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It is known how to map formally the noncommutative theory into a matrix theory 
by the Weyl map. Let x M and y v be Hermitean operators satisfying [x^,y u ] = iO^, 
which is the operator analogy of (4). One can define the Weyl map by 



= y d 2 x(j){x) A(x), A = J _^A e ^"e- to . 
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and in this way the action (1) can be represented (formally, with a suitable normal- 
ization of the trace and suitable assumptions concerning the fall off of the functions 
4>(x) at infinity) as 
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where is an anti-hermitian derivative which together with x^ satisfy 

The operators (or matrices) in (7) are all infinite dimensional. However, the as- 
sociated Heisenberg algebras of exponentials of x^ (and <9 M ) allow finite dimensional 
representations, as already noted by Weyl. This is the key to formulating the non- 
commutative theory on a finite periodic N x N lattice with lattice spacing a (see 
[5] for a more general formulation). The algebra (7) is replaced by 

Z„Z V = e- 2m ^Z v Z^ D^D^ = e 2 * iS ^ N Z„, (8) 

where 

f- ' f 1 ' w— a 2 N 2 ' 

Periodicity on the lattice fixes the quantization of momentum and 9 

27rm u n na 2 N ( ^ 2n\ 

where G 0, . . . , N — 1, and n is an integer. In the following we choose n — 1 for 
simplicity. 

A simple realization of the above scenario is obtained by choosing Z M , fi = 1, 2 

to be the standard N x N shift and clock matrices {Z\)jk = Sj+i,k, {^2)jk — 
( e 4*i/Ny-i Sjk ^ f )i = ^t)(jv+i)/ 2j ^2 = (^ 1 )(A f +i)/2 j the lattice size TV being odd. 

In this way the lattice version of (3) becomes 

= ^ X>(y)&(*) ^w-'K (ii) 

and the operator A becomes an N x N matrix 

N 

A(x) = (z^Z™ 2 e~™ 0mim2 ) e ~ 2 ™ m ^/ aN . (12) 

mi ,rri2=l 

It generates an explicit map from the lattice field <j)(x) defined at the N x N lattice 
points x to that N x N matrix defined by (5), 

^^E^O^*)' #*) = -^M(z) (13) 
such that the action is given by (6) with N x N matrices and the partition function 



is 



Z(R,G) = J d0e~ s M, (14) 



S[(j)} = iVTr 
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(15) 



The dimensionless lattice parameters are 



R = l -m 2 a\ G = l -g 2 a 2 . (16) 

The ordinary (commutative) 4 theory on the lattice has two different continuum 
limits when the space-time volume goes to infinity, one "trivial" associated with the 
fine tuning the approach to R = and G — 0. It corresponds to the standard c/> 4 
two-dimensional field theory where only normal ordering is needed to renormalize 
the theory. The other limit is associated with the phase transition between a phase 
in which the symmetry <fi — > — </> is broken and a symmetric phase. It occurs for 
finite negative values of R and finite positive values of G and the continuum limit 
can be identified with the continuum limit of the Ising model, i.e. with a c— 1/2 
conformal field theory. 

In the noncommutative case the continuum limit of course requires an infinite 
lattice, i.e. we have to take N to infinity. In addition we are interested in a limit 
in which the lattice spacing can be viewed as scaling to zero while the dimensionful 
parameter 9 stays finite. Two possible scenarios appear possible: 

If non-commutativity can be considered as a perturbation imposed on a commu- 
tative theory then a kind of "double scaling limit" of the matrix model (15) suggests 
itself: 

tt9 

N-^oo, a 2 = —, a 2 N 2 = 7i9N. (17) 

Thus the physical space-time diverges like N 1 and we can write the matrix model 
action (15) as 



S[<j>\ = NTr 
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(18) 



The noncommutative parameter 9 only appears in the dimensionless combinations 
m 2 9 and g 2 9 and implicitly in factors of iV appearing in (18). This type of continuum 
limit seems to be most naturally associated with the trivial gaussian fixed point of 
the commutative theory. Note that this "double scaling" limit iV — > oo is quite 
different from the usual planar limit of matrix models and thus also quite different 
from the double scaling limit of the 2d matrix models studied in the context of 2d 
quantum gravity. 

If, on the other hand, non-commutativity alters the commutative theory pro- 
foundly one should follow the general strategy of non-perturbative lattice stud- 
ies in order to identify a continuum limit: identify suitable physical observables 
usually related to a divergent correlation length of the lattice field. Fine-tuning 
the coupling constants such that the observables are constant when expressed in 
terms of the "physical" correlation length a ■ £, where £ is the correlation length 

1 For other limits which allow for a finite physical space-time see [7, 6, 5]. 
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in lattice units, will then lead to a relation between the coupling constants and 
the lattice spacing a eg. a ~ (R — R c ) v and allow us to identify the (potentially) 
non-perturbative renormalization of the coupling constants required to obtain a con- 
tinuum theory. We can then potentially investigate another "double scaling" limit 
in which N (R — R c ) 2u ~ 6 is held constant as iV — > oo and R — > R c . Such a phase 
transition would generalize the usual Ising transition of the commutative theory to 
the noncommutative case. 

The first scenario has the attractive feature that we know ahead of time the value 
of the non-commutativity parameter 9 and the dependence of lattice spacing on the 
matrix size N. We will show numerical evidence which supports scaling according 
to this scenario although we stress that the existence of such a double scaling limit 
does not rule out the possible existence of additional continuum limits based on the 
second scenario. 

Susceptibilities (i.e. correlation functions integrated over all space-time points) 
have been used in conventional lattice field theory to identify the couplings asso- 
ciated with a divergent correlation length, and the critical indices associated with 
the susceptibilities lead to a classification of the continuum Euclidean field theories. 
Susceptibilities are suitable observables also in the noncommutative case. In fact 
one has 

J d d x 4>{x) * <f)(x) = J d d x <j>{x)<j>{x), (19) 

so they are maximally close to the commutative observables. Yet they may contain 
additional interesting information about the infrared behavior of the noncommuta- 
tive propagators. In the computer simulations to be reported below we will indeed 
use (generalized) susceptibilities to identify potential phase transitions. 



(P$(-P)> 



Physically, they are just the (mod square) of the Fourier component of the noncom- 
mutative field. 



Stripes 

The possibility of having a ground state which breaks translational invariance was 
first observed in [1]. The quadratic part of the effective action, obtained from (1) 
by a one-loop self-consistent Hartree calculation, is 

rW(p)=p* + m l + A 9 HnAl, K= e ^ l +l/ ^ ^ 

A denotes a UV cut-off, m 2 R the renormalized mass term and A is some constant. 
The mixing of IR and UV modes is manifest in the form of A p . In addition the 
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momentum dependence of A p implies that r^ 2 ^(p) can have a minimum for 



Vl = Ag 2 - 



1 



(21) 
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For sufficiently large coupling this analysis would predict that modes with momen- 
tum p = p c would be the first to become unstable as the mass parameter uir were 
tuned to zero. Thus any phase transition would then be driven by condensation of 
a non-trivial momentum component of the field. This is to be contrasted with the 
usual Ising transition in which the p = component of the field condenses to yield 
the broken symmetry phase. Using (16) and (17) this result can be written in lattice 
variables 



Note that for small G the matrix size N has to be large (N > 2ir/\/A G) in order 
for (22) to have a solution. Conversely, we might expect an exotic vacuum state to 
exist for any coupling in the continuum limit N — > oo. 

In two dimensions it is usually difficult to have a spontaneous breaking of a 
continuous symmetry like translation invariance. The fluctuations associated with 
the corresponding massless Goldstone boson will spoil any long range order. How- 
ever, the arguments leading to this conclusion are strictly speaking not valid for 
non-local field theories and there exist explicit examples of non-local field theories 
in two dimensions which exhibit spontaneous symmetry breaking 2 Since noncommu- 
tative field theories are non-local and show a mixing of IR and UV properties one 
cannot rule out that the simple one-loop arguments given above in favor of a non- 
translational invariant ground state could still be valid. In fact the renormalization 
group treatment in [2] predicts this ground state even in two dimensions. 

In the following we will look for signals indicating the appearance of a non-trivial 
ground state in the lattice system (15). 

The computer simulation 

In these first explorative simulations we have used a simple metropolis algorithm for 
stimulating the matrix system (14)-(15). Specifically, the matrix field is updated by 
the addition of a random hermitian noise matrix rj 



We tune e with N so as to achieve an acceptance rate close to 50%. Most of our runs 
employed 5 x 10 6 such matrix updates for a given set of values of the parameters m 
and g. If we only consider the quadratic part of (15) the system is equivalent to an 
ordinary gaussian lattice field theory via the map (13) and this has been a useful 
check of the computer code. 

2 One example is the so-called crumpling phase transition in crystalline membranes. 




(22) 
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Let us now turn to the actual simulation with the noncommutative 0* 4 term. We 
have looked in the (m,g) coupling constant space for signals of a nontrivial vacuum. 
We have not attempted to map the entire phase diagram but concentrated instead 
on the slice 9 = 1.0 and g 2 /4 = 1.0 with varying physical mass parameter to. In Fig. 
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Figure 1: the average of the (mod square) of the Fourier amplitude for 4 lattice 
modes vs mass parameter m 2 (for lattice size N — 9). The 4 modes correspond to 
momentum states (24) where (to, n) is (0,0), (1,0), (1, 1) and (2,0). 

(1) we show the generalized susceptibility Xk corresponding to the (mod squared) 
Fourier component of the field 0(x) on an N = 9 lattice for various values of the 
lattice momentum 

( 27im 2nn\ . , . . . . 

(*i,*2) = (^r>^rj> (p = */«)■ (24) 

One sees a clear condensation of the usual p = (0, 0) mode for values of to 2 < —6. 
Contrast this with the Fig. (2) which shows the same set of (unsubtracted) suscep- 
tibilities for iV = 23. Now for to 2 < —6 we see that the dominant susceptibility 
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Figure 2: the average of the (mod square) of the Fourier amplitude for 4 lattice 
modes vs mass parameter m 2 (for lattice size N = 23). The 4 modes correspond to 
momentum states (24) where (m,n) is (0,0), (1,0), (1, 1) and (2,0). 

corresponds to the lattice (1,0) mode - the vacuum in this region of the parame- 
ter space shows a spatially varying order parameter as predicted by the one loop 
calculation - a striped phase. Since the system still possesses a symmetry under 
exchange of axes we find that the condensation takes place along either the 1-axis 
or 2-axis randomly as the mass is varied. Thus, in our plots we show the average 
< \4>(n, m)\ 2 + \<f>(m, n)\ 2 >. Another way to illustrate graphically the existence 
of this phase is by considering the transverse-averaged correlator in position space. 
Fig. (3) shows a plot of < 0(xi)0(O) > for m 2 = —10.0. Notice that the correlator 
shows one complete oscillation within a lattice length. Since we are deep within 
the striped phase this correlator will essentially measure the behavior of the order 
parameter directly. Along the second direction the correlator is positive definite as 
shown in Fig. (4). This behavior is what we would expect of a striped phase. For 
the smaller lattice A = 9 we do not see any asymmetric behavior in the form of this 
position space correlator in the two directions, and indeed, as Fig. (5) indicates, the 
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Figure 3: (0(xi)0(O)) vs. x x for m 2 = -10.0 and iV = 23 

correlator is of conventional (positive definite) form over the entire lattice. This is 
consistent with the idea that for small N, there are no solutions to (20) and the 
minimum of the quadratic effective action is still found for mode (0, 0). 

We have also examined the system for N = 39. Fig. (6) and Fig. (7) show the 
corresponding correlators along the two lattice axes for m 2 = —10.0. From these 
two pictures we infer that the condensed mode is now (2,1). Thus, rather than 
stripes we actually see a checkerboard-like arrangement for the vacuum state in this 
case. The fact that the lattice mode number increases with larger lattice size for 
fixed physical parameters is encouraging as this is the behavior that is required if 
this effect is to survive the large volume limit N — > oo. Indeed, (20) indicates that 
the mode number n should increase like n ~ for N — > oo to ensure that the 
physical momentum at which condensation occurs remains fixed. 

Discussion 

We have shown that it is indeed possible to observe and investigate non-trivial 
phenomena like the appearance of "stripes" in noncommutative field theory within 
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• • m 2 =-10.0,g 2 /4=1.0 
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Figure 4: (0(x 2 )0(O)) vs. x 2 for m 2 = -10.0 and N = 23 



the non-perturbative framework of lattice simulations. 

We have here looked at the very simplest model, a two-dimensional scalar field 
theory simply to verify that phenomena like the appearance of a "striped phase" can 
be observed. Our data show the appearance of such a phase in an appropriate double 
scaling limit. Notice that the issue of possible I.R divergences in two dimensions does 
not preclude such a study - the lattice regulates any such divergences automatically 
and the issue returns only when the continuum limit is taken. Our current data do 
not address this issue - much larger volumes must be simulated in order to determine 
whether this phase truly survives the large N limit. Nevertheless, we feel that our 
simulations lend nonperturbative, albeit qualitative support to the analysis of [1] 
and [2]. However, it is clearly of more interest to do this in higher dimensional 
theories. The lattice formalism has been developed [5] and just awaits use. Indeed, 
while this paper was in preparation we received a preprint [8] in which the same 
model was examined in 2 + 1 dimensions with similar conclusions. 
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